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Abstract 

Order /3qO^ +1 QCD corrections to top quark decays into W + and H + 
bosons are computed to all orders in perturbation theory. Predictions for 
the radiative corrections to the top quark width are compared with the es- 
timates from BLM scale setting procedures. The results of the summation 
are shown to greatly improve understanding of higher order corrections in the 
limit m\Y, mn — ► rat, where the BLM scale setting method is known to fail. 
Attempts to reduce nonperturbative error by substituting the running mass 
for the pole mass in the expression for the decay are shown to fail in the limit 
mw, rriH —* m% because of subtleties in the treatment of phase space. 
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Now that the top quark has been observed by CDF |]J and DO 0, future experiments will 
attempt to do precise measurements of its properties. To date, only the leading order QCD 
corrections to the decays of the top quark into W + and H + bosons have been computed ||(4] . 
The theoretical predictions suffer from considerable uncertainty since the renormalization 
scale for the strong coupling is undetermined at leading order. Corrections of 0{(3oa 2 s ) are 
known |||6| and can be used to set the scale in the leading order calculation via the scale 
setting method of Brodsky, Lepage and MacKenzie(BLM) H. (/3q is the first coefficient of 
the QCD beta function: (3q = (11 — 2/3 rif)/Aix.) In a previous letter 0, I showed that the 
BLM scale setting method, when applied to the decays t — > W + + b and t — > H + + b, gives 
unreasonably small scales in the limit mw, mn m t- Here I will perform a summation 
of 0(/3o«" +1 ) terms in the perturbation series and find that the summation gives a much 
more reasonable picture of the effect of higher order corrections when the BLM scale setting 
method fails. For the case of t — > W + + b, m w <C m t , so the summation of 0{(3qO^ +1 ) 
terms affords little improvement over BLM scale setting. However, the all order summation 
of 0(/3q«™ +1 ) terms could be phenomenologically relevant if it turns out that the top quark 
can decay into a new heavy particle such as a charged Higgs. More generally, the BLM scale 
setting method is widely used as a tool for estimating the impact of higher order corrections 
on perturbative QCD predictions. For this reason, it is important to understand when it 
can fail and what means are available for improving estimation of higher order corrections 
when it does fail. 

In the BLM scheme, the scale is chosen so as to absorb the 0{f3oa 2 s ) term into the leading 
order correction: 

Tit - W + /H + + b) = T (l + (r Q + ( ri/3o + S)a s + ...)) 

= T (l+ 2a ^ BLM \ r + 6a s + ...)) (1) 



3tt 



where 

[J>blm = mt exp 



-n 

2r n 



r (^^ + 6) = W ^fi + 24)fi-4) , 

8V27T V m t ) \ m l t j 

r (t ->H+ + b) = \v tb \ 2 ^^cot 2 p fi-^fV. 

8a/27t V m t J 

(In these expressions, as well as throughout this paper, I will neglect the mass of the b 
quark.) The resulting BLM scale turns out to be extremely small, especially as and 
mil approach m t . For example, in the case of a 175 GeV top quark decaying into a 150 
GeV Higgs and massless b quark, the computed BLM scale is 200 MeV. Taken at face value, 
the BLM scale tells us that the radiative corrections to this process are uncontrolled even 
though the energy released in the decay, 25 GeV, is much greater than the QCD scale. It 
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is not clear that the extremely small BLM scales reflect a genuine physical effect, since the 
very large scales that result from BLM scale fixing appear to be a consequence of a small 
leading coefficient in the perturbative series. 

Since the scale dependence of the QCD prediction is an artifact of the finite order trun- 
cation of the perturbative series, any scale fixing scheme is essentially a guess of the size 
of higher order coefficients in the perturbative expansion. Inserting the one-loop expression 
for the QCD coupling constant, 



into Eq. ([]]) and expanding in powers of 0(/3Qa s (m t ) n+1 ), one finds that the series after 
BLM scale setting is equivalent to: 



The BLM scale setting scheme can be thought of as a crude summation of 0(/9o«s +1 ) terms 
in the perturbative series with a simple ansatz for the coefficients of these terms. This scale 
setting scheme suffers from obvious deficiencies. First of all, the size of these coefficients can 
be grossly overestimated (or underestimated) if the leading (or next-to-leading) coefficient 
is anomalously small. This is what appears to be going on in the case of top decay. In 
Ref. || it is observed that the rapid decrease of the BLM scale as m w , m H — > m t appears to 
be a consequence of the leading order term becoming small rather than the next-to-leading 
term becoming large. It is not obvious that the small leading coefficient implies the kind of 
growth of higher coefficients predicted by the BLM scale setting scheme. Computation of 
higher order 0(/?q«" +1 ) terms allows one to assess the perturbative series as a whole so as 
to be less sensitive to the vagaries of a single coefficient. The second deficiency concerns the 
asymptotic behavior of the series. The 0(/?q«™ +1 ) term of the perturbative series contains 
a piece proportional to n^a™ +1 which comes from the insertion of n fermion bubbles in the 
gluon propagator of the leading order correction. Direct computation of these graphs shows 
that their coefficients exhibit factorial growth at large orders, whereas the BLM scale setting 
scheme predicts that the coefficients of higher order terms in the series grow geometrically. 
This geometric series diverges if tx/tq = l//5oo; s . This divergence is an artifact of the BLM 
scale setting method and the real 0(/?Qa" +1 ) series does not exhibit any such divergence. Of 
course, the true 0(/?Qa™ +1 ) series is asymptotic and so one must be careful to define what is 
meant by its sum. This is done by means of a Borel transform, as will be discussed below. 
The result of this summation contains no spurious divergences. 

It is important to acknowledge some of the shortcomings of the summation of 0(/?o«s +1 ) 
terms in the perturbative series. There is no physical limit of QCD in which these terms 
dominate. Nevertheless, one may expect these terms to make up a large part of the 0(a™ +1 ) 
correction because of the factor (11 — f ?i-/) n . Indeed, many next-to-leading order corrections 
in perturbative QCD calculations are dominated by the O^qo 2 ) term, and so the compu- 
tation of this piece gives a decent approximation for the complete 0(a 2 ) correction. (For 



ol s \IJ>blm) 



a s (m t ) 



1 + f3 a s (m t )\\i(^ 2 BLM /m 2 t ) ' 
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examples of this, see Refs. JF]||.) At higher orders, the coefficients of the 0{(3qO^ +1 ) terms 
grow like 2 n n\. The numerical enhancement of these terms suggests they may be a useful 
guide for estimating the size of higher order corrections. 

The decay rates for t — > H + + b and t — > W + + b can written in the form: 

F (t W +/ H + + h) = r (\ + 2 _^A g a s {m t nr n (3 n Q + S^j (3) 

The S n are 0(/3q or smaller and will not be computed. The series J2^=o r n (Po^ s ('mt)) n 

does not converge because r n ~ 2 n n\ for large n. The sum of the infinite series is defined by 
its integral representation in terms of its Borel transform: 

1 rco 

S^^— due- u ^B(u); (4) 
fj a s Jo 



n=0 "- 

The integral representation for Soo is still not well defined because the function B(u) has 
poles on the positive real u axis. These are the infrared renormalons. In order to define the 
integral, the contour must be deformed into the complex plane. The result of the integral 
has an (imaginary) ambiguity depending on whether one chooses to deform the contour 
above or below the real axis. If the Borel transform of the perturbative series has a pole 
at u — Uq, then the ambiguity in the integral representation for the series, which is the 
difference between going above or below the real axis, is just the residue of the pole at u$: 

Im ~ e -" o//3oQs(mt) ~ ( ^£ ) . (5) 

\ m t J 

The ambiguity in the sum of the series has the same size as a power correction one would 
expect to receive from non-perturbative QCD. It is customary to define the sum of the 
perturbative series as the principal value of the Borel integral and use the imaginary part of 
the Borel integral as an estimate of the uncertainty in the perturbative theoretical prediction. 

It is important to understand that the error estimated this way is usually an underes- 
timate of the true size of nonperturbative effects. One way to see why this is the case is 
to note that Aqcd appearing Eq. (|5|) is the parameter in the one- loop expression for the 
running coupling, which is ~ 100 MeV, while the intrinsic scale of nonperturbative QCD 
effects is typically closer to 1 GeV. For observables with an operator product expansion, 
the size of the uncertainty in the perturbative series as defined by the imaginary part of the 
Borel transform usually underestimates the size of the lowest dimension hadronic matrix 
element which makes up the leading non-perturbative correction 0. 

For processes in which there is no gluon self coupling at leading order, the Borel transform 
can be compactly expressed in terms of the one loop correction computed with finite gluon 
mass-squared, A 2 |10| : 

7vu Jo \mf J 
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where r' (X 2 ) = dro(X 2 )/dX 2 . In this expression, C is a renormalization scheme dependent 
constant which turns out to be —5/3 in the MS scheme. The coefficients r n can be computed 
by taking derivatves of the Borel transform: 



d n 



(7) 



Inserting the expression for the Borel transform in Eq. (^j) into Eq. (H) results in the following 
expression for 



00 ^(A 2 K(A 2 ) + ro(Ai i- ro(0) . 
(I p a s 



In this expression, 



(8) 



$(A 2 



1 



f3 a s V 7r 



-tan 1 



1 + (3 a s \n(X 2 e c /m 2 ) 



and 



Ai 



-m t exp 



/3 a s (m t ) 



C 



e~ c A 2 



QCD- 



The imaginary piece of the Borel transform comes from the term proportional to tq(X 2 l ). 
The uncertainty in the summation of the series is defined to be: 



5S n 



Imr (A|) 
7ip a s 



(9) 



For the decays t — > W + + b and t — > H + + 6, the principal value of the Borel integral is 
approximately the same as the partial sum of the series truncated at its smallest term, and 
the error given by Eq. (|^) has nearly the same numerical value as the smallest term in the 
series. 

In Fig. [I], the result of applying the BLM scale setting method to t — > W + + b is compared 
with the all orders summation of 0(/?q«" +1 ) terms. I define T = r + ST, and plot ST/T 
as a function of rriyy/m 2 . After applying the BLM procedure, ST/T is given by: 



6T 

r~ 



BLM) 



37T 



2a s (m t ) 



r 



3tt 1- Aja 8 (mt)(ri/ro)' 



(10) 



After summing all 0(/?q«™ +1 ) terms in the perturbation series, <5r/r is instead: 



<5T _ 2a s (m t ) 

r 3tt 



s a 



(ii) 



The value of fl a s (m t ) is taken to be 0.066. For m^/m 2 < 0.7, the BLM scale setting 
method gives a result that agrees fairly well with the all orders summation of 0(/3Qa™ +1 ) 
terms. For larger values of trw, the BLM scale setting method ceases to give meaningful 
results because Eq. fllOD has a pole near m^/m 2 = 0.85. As discussed earlier, this spurious 
pole is an artifact of the geometric approximation of the full 0(/?o a ™ +1 ) series. When the 



5 



higher order terms in the series are computed correctly the sum of the series exhibits no 
pole. Perturbative corrections never exceed 17% for all values of my/. 

The uncertainty of the summation due to the asymptotic nature of the series is studied 
in Fig. ||, where I plot SSoo/S^. As expected, the uncertainty increases as my/ — > m t , but 
the uncertainty in the calculation of Soo is never greater than 20%. Since this enters the 
decay rate multiplied by a factor of 2a s (mt)/37T = 0.023, this amounts to an uncertainty 
in the prediction for the total rate that is less than 0.1%. As mentioned earlier, true non- 
perturbative corrections are expected to be larger. Nevertheless, it is safe to conclude that 
perturbation theory works well for top quark decay all the way to the kinematic limit. This 
conclusion is in stark contrast to what I find based on the application of the BLM scale. 
For m^/m 2 = 0.8, \iblm — 484 MeV. One could be fooled in to thinking that perturbative 
QCD does not work. Examination of higher order 0(/?q«™ +1 ) terms indicates that it does. 

For the physically relevant case of m^/m 2 = 0.21, I find 



The estimate for the error is Aqcd /i?k with Aqcd = 1 GeV. This gives a significantly higher 
error estimate than what results if SSoo is used. The reasons for this more conservative error 
estimate were given earlier. 

The results for the decay t — > H + + b are similar to the results for t — > W + + b decay. In 
Fig. H and Fig. |], I plot <5r/r and SS^/Sco. The BLM scale setting method breaks down 
for m 2 H /m 2 > 0.6, as opposed to my/jm\ > 0.7 in the case of W decay. The total correction 
to the decay rate is roughly 15% for all values of the Higgs mass. Uncertainties associated 
with the summation are approximately the same size as in the case of t — > W + + b, always 
less than 0.1% of the total decay rate. Again, perturbation theory works for all values of 
m#, despite the failure of the BLM procedure for large m#. 

The piece of the next-to-leading order corrections that has not been computed is expected 
to be 0(a s (mt) 2 /tt 2 ) ~ 0.001. This is actually smaller than the leading nonperturbative 
corrections of 0(Aqc , d/"t^)- Unless the uncertainty due to nonperturbative corrections 
can be eliminated, computation of higher order corrections to top decay will afford no 
improvement of the theoretical prediction. 

In Ref. 0, it is argued that rewriting the perturbative series in Eq. (|XJ) in terms of the 
top quark running mass will reduce the sensitivity of the perturbative series to long distance 
scales. The series relating the top quark pole mass to the top quark running mass (in the 
MS scheme) is given by 



where d = 4, d\ = 71/8 + n 2 , and higher order di can be found in Ref. [M. The Borel 
transform of this series suffers from a renormalon pole at u — 1/2, just like the series for 
top quark decay. When the top quark pole mass is traded for the running quark mass, the 
poles at u = 1/2 cancel and the leading pole in the Borel plane is located at u = 3/2. 
Consequently, the error in summation of the series is 0(Aq CD ). 

Replacing m v t ole with the series invloving m7(m 4 ) results in the following expressions for 
top decay: 



6T 

r~ 



-0.135 ± 0.006. 




(12) 
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F(t ^W+ + b)~ m 3 (l + 2w) (1 - wf (l + £(/3o«,) n (r™ + n T + o W ^ 2W2 \ d n) + 



n=0 



2(l + 2w)(l -w) 



F(t -> F+ + 6) ~ mf (1 - />) 2 (l + (r? + J^U n ) + , (13) 



n=0 



2(1 — /i) 



where u> = m w jra 2 and = m 2 H /mi 2 . Note that in Eq. (14) of Ref. 0, the expression for 
the series after trading the top quark pole mass for the top quark running mass is incorrect. 
Consequently, the calculations of the BLM scales for the series involving the running top 
quark mass in that paper are wrong. When the correct expressions in Eq. are used, the 
BLM scales are roughly 0.1 — 0.2 m t and exhibit only a weak dependence on mw, mu- 

After substituting the pole mass for the MS mass, I find the 5T /To shown in Fig. |5|. In 
this case, the BLM scale setting method gives almost identical results as the all-order sum- 
mation of 0(/3o<a;" +1 ) terms. The uncertainty in the summation 5Soc is negligible compared 
to uncomputed 0(a 2 ) corrections for all values of mw, mu- The sum is dominated by the 
leading order coefficient in the perturbative series, the next-to-leading order term gives a 
small correction, and higher orders are negligble. This is why the BLM setting scheme gives 
such a good approximation to the summed series. 

For mw <C m t , the prediction for the radiative correction is roughly the same whether one 
uses Eq. ( O) or Eq. (|3|). As mw — ► mi, the prediction for the radiative corrections becomes 
much larger if the series is expressed in terms of the running mass. This is because the phase 
space for the decay is incorrectly given in terms the running mass. The decay rate given by 
Eq. (|l^) vanishes when mw, mu = m~i instead of mw, mu = m\ ole . Trading the pole mass 
for the running mass results in an unphysical pole at mw, mu = mi in the expression for 
the radiative correction. Therefore, despite the formal reduction of non-perturbative errors, 
this substitution is not reasonable near the kinematic limit. 

Even when mw <C m t , it is questionable whether rewriting the decay rate in terms of the 
top quark running mass will really improve theoretical predictions for the top quark width. 
It is certainly true that the perturbative series for radiative corrections to the tree level 
rate is improved, in that terms higher order than n = 2 in 0(/3Qa™ +1 ) series are rendered 
negligible. The price to pay for this is that the tree level rate is now expressed in terms 
of a variable whose relation to physical observables suffers an O(Aqcd) ambiguity. The 
top quark mass experimentally defined as the peak in the invariant mass distribution of 
a W boson and b jet in a high energy hadronic collision is clearly the pole mass, up to 
corrections of 0(Aqcd)- If one tries to use Eq. (fL2) to determine the running quark mass 



from the experimentally determined top quark mass, the value of the running top quark 
mass will be subject to O(Aqcd) corrections, both from the fact the experimental top quark 
mass is not the pole mass, and because the series relating m? t ole and m ] f s {m^) suffers from 
an O(Aqcd) ambiguity. Hence, the prediction for the top quark width will still recieve 
O(Aqcd) corrections. Only if several observables are expressed in terms of the running 
mass, then the running mass eliminated to derive relations between these observables, will 
there be any hope of deriving theoretical predictions free from O(Aqcd) uncertainties. 

In summary, I have discussed the improvement in the estimation of higher order correc- 
tions to the top quark decays t — > W + + b and t — ► H + + b. The BLM scale setting gives 
meaningless results as mw, mu m t , but summation of 0(/5ga™ +1 ) corrections shows that 
perturbation theory is well-behaved even in this limit. 
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FIG. 1. ST/T for t -> Solid line - results after BLM scale setting; Dashed line - 

summation of 0(/3qQ!™ +1 ) terms 
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FIG. 3. 5T/Tq for t — > Solid line - results after BLM scale setting; Dashed line - summa- 

tion of 0(/^a™ +1 ) terms 
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FIG. 5. 5T/Tq after top quark pole mass is replaced with the running mass (MS scheme). Solid 
line- t -> W + + b; Dashed line- t^H + + b line. 
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